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ABSTRACT 

We report the discovery of a debris system associated with the ~ 30 Myr 
old G3/5V star HD 12039 using Spitzer Space Telescope observations from 3.6 
- 160/rni. An observed infrared excess (Lir/L* = 1 x 10~ 4 ) above the expected 
photosphere for A > 14/xm is fit by thermally emitting material with a color 
temperature of T~ 110 K, warmer than the majority of debris disks identified 
to date around Sun-like stars. The object is not detected at 70/im with a 3cr 
upper limit 6 times the expected photospheric flux. The spectrum of the infrared 
excess can be explained by warm, optically thin material comprised of blackbody- 
like grains of size > 7/mi that reside in a belt orbiting the star at 4-6 AU. An 
alternate model dominated by smaller grains, near the blow-out size a ~ 0.5/im, 
located at 30-40AU is also possible, but requires the dust to have been produced 
recently since such small grains will be expelled from the system by radiation 
pressure in ~ few xl0 2 yrs. 
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1. Introduction 

Dust in the inner Solar System can be seen via scattered light as the zodiacal glow visible 
to the naked eye near sunrise and sunset in dark clear skies. This dust is located within 
the terrestrial planet zone, arising primarily from collisions among rocky debris (asteroids) 
between Mars and Jupiter and secondarily from comet ejecta. Dust produced by collisions 
between remnant planetesimals is also expected to exist within the Kuiper Belt beyond 
Neptune. The region between 5 and 30 AU in our solar system is kept relatively free of 
dust due to the dynamical action of the giant planets (Moro-Martm & Malhotra 2002, 2003; 
Moro-Martin, Wolf & Malhotra 2005). Comparable planetary systems around other stars 
might be expected to produce similar dusty debris belts and gaps. 

Examples of thermal emission from debris disks around fairly luminous main-sequence 
stars, including the proto-type Vega (Aumann et al. 1984), were identified with IRAS (e.g., 
Walker & Wolstencroft 1988; Backman & Paresce 1993; Mannings k Barlow 1998) and ISO 
(e.g., Habing et al. 2001; Spangler et al. 2001; Laureijs et al. 2002). Most of these systems 
are associated with A-type stars since neither observatory had the sensitivity to detect the 
amount of material associated with Vega around low-luminosity Sun-like stars at distances 
beyond about 15-20 parsecs. Most debris systems identified to date are characterized by cool 
(T ~ 50 — 90K) disks, suggestive of material at gas giant or Kuiper Belt distances from their 
parent stars, r^ust ~ 10 — 100 AU (see, e.g., Backman & Paresce 1993; Decin et al. 2000, 
2003; Lagrange, Backman & Artymowicz 2000; Zuckerman 2001 and references therein). 

Since much of the zodiacal dust in our solar system arises from collisions in the asteroid 
belt (e.g., Grogan et al. 2001; Nesvorny et al. 2002, 2003), we might expect to observe excess 
emission from debris in a similar zone around other stars (Stern 1994; Zhang & Sigurdsson 
2003). However, only a handful of debris disks associated with less luminous, Sun-like stars 
have been found to exhibit peak infrared excess emission from warm dust (T ~ 100 — 300K) 
suggestive of terrestrial zone material (r dust ~ 2 — 5 AU). Herein, we report the discovery 
of another example that is associated with HD 12039 (DK Cet, HIP 9141), a young («30 
Myr) G3/5V star (Houk et al. 1988) at 42±2pc (Perryman 1997). The infrared excess 
was identified in the initial enhanced data product release from our Spitzer Space Telescope 
(Spitzer. Werner et al. 2004) legacy science program titled FEPS (Formation and Evolution 
of Planetary Systems: Meyer et al. 2004, 2005; Hines et al. 2004a). 
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Sections 2 and 3 briefly describe the FEPS program and the Spitzer observations of 
HD 12039. Section 4 describes our discovery of the debris system associated with HD 12039, 
and presents some of the properties of the star and its infrared excess. We further characterize 
and model the debris system in Section 5. Section 6 explores the origin of the emitting 
material and how HD 12039 fits into models of disk evolution and planet formation. We also 
compare HD 12039 to other warm debris systems associated with Sun-like stars. 

2. The FEPS Program and VI. 1 Sample 

The FEPS Spitzer Legacy program encompasses observations of 328 Sun-like (~ 0.8-1.5 
M ) stars with ages ranging from ~ 3 Myr to 3 Gyr. The survey enables us to examine 
the frequency of stars that have circumstellar gas and dust, as well as to infer some of the 
properties of detected dust (Meyer et al. 2004; Meyer et al. 2005). Our ultimate goal is 
to assess whether debris systems such as that found in our solar system are common or 
rare around Sun-like stars in the disk of the Milky Way. FEPS uses all three Spitzer science 
instruments to provide spectral coverage from 3.6/xm to 70/zm, and includes A ~ 7— 35;um low 
resolution IRS spectra. We also obtain 160 fim observations for a subset of the FEPS stars 
located in areas of low infrared cirrus background. In addition, FEPS provides models of the 
photospheric emission that have been fit to available optical and near-infrared photometry 
for each star. 

The initial FEPS data release contains ~ 10% (33 stars) of our total sample. A de- 
tailed description of the observing strategy, data analysis and best-fit, Kurucz photospheric 
emission models for these stars are provided in the FEPS Explanatory Supplement (Version 
1.1: Hines et al. 2004a) 1 . The Spitzer data and documentation for the FEPS program are 
available directly from the Spitzer Legacy Science Archive 2 . 

The 33 stars were randomly chosen by the SSC for scheduling between UT May 2004 
and UT September 2004. They cover a wide range of properties within the selection criteria 
used for the entire FEPS program, and therefore provide a good sub-sample from which 
to begin looking for interesting objects. In particular, these stars cover a broad range of 
effective temperatures 5121K < T 33 < Q227K (compared with A299K < T FEPS < Q769K 



1 Subsequent to the release of VI. 1, the IRS spectra were reprocessed by an updated SSC pipeline that 
significantly improved the calibration. The MIPS Instrument Team also released an improved data analysis 
pipeline that reduced the noise at 70/zm. These updated IRS and MIPS data are used herein and have been 
released publicly in FEPS V2.0. 

2 http://data. spitzer. caltech.edu/popular/feps/ 
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for the entire sample), have an age range 7.2 < log(age) < 9.2 [yrs], and lie at distances 
20 < d < 127 pc. FEPS stars with younger ages 6.2 < log(age) < 7.2 and greater distances 
d > 130 pc are primarily cluster members and are less well represented by these 33 stars. 

FEPS results to date include: 1) the frequency of young, optically thick disks detected 
with the IRAC instrument (Silverstone et al. 2005); 2) cool debris systems having excesses 
detected only at MIPS 70/im and therefore comparable to our own Kuiper Belt (Meyer et 
al. 2004; Kim et al. 2005): and 3) constraints on the amount of remnant gas in the disk 
surrounding HD 105 (Hollenbach et al. 2005). We have also investigated the occurrence of 
debris disks in the 100 Myr old Pleiades star cluster (Stauffer et al. 2005). We focus the 
current study on the Sun-like star HD 12039, which exhibits dust emission with temperatures 
that are warmer than found in typical disks around other stars. 

3. Spitzer Observations of HD 12039 

In this section we briefly review the observing strategy and data reduction methods 
specifically for HD 12039. Observations of the other 32 stars were executed and processed 
in a similar manner. Additional details are presented in the Explanatory Supplements for 
the FEPS VI. 1 and V2.0 data releases (Hines et al. 2004a, 2005), and in Kim et al. (2005), 
Bouwman et al. (2005), and Silverstone et al (2005) for MIPS, IRS and IRAC respectively. 

Multi-band Imaging Photometer for Spitzer (MIPS: Rieke et al. 2004) observations of 
HD 12039 were obtained using the small field photometry mode on UT 2004 July 11 at 24, 
70 and 160/im. Two cycles of 3s Data Collection Events (DCEs) at 24/im, two cycles of 10s 
DCEs at 70/xm, and four cycles of 10s DCEs at 160/xm were obtained. After initial processing 
by the SSC S10.5.0 pipeline to provide reconstructed pointing information, the MIPS data 
were further processed using the MIPS Data Analysis Tool (DAT, ver. 2.9; Gordon et al. 
2004, 2005). This includes the "enhancer" routines that correct for distortion in individual 
images and combine the images onto a sub-sampled mosaic. Aperture photometry using 
IDP3 (version 2.9: Schneider & Stobie 2002) was performed on the individual 24/im DCE 
images, and on the mosaic images at 70/im and 160/im. We adopt 14'/7, 29'/7 and 48'/0 target 
aperture radii at 24, 70, and 160/im, respectively. Background annuli from 29 / .'4 - 41"7, 39'/6 
- 79'.'2 and 48'.'0 - 8(f!0 were used for the three bands respectively. The median background 
per pixel was scaled to the appropriate target aperture size and subtracted from the summed 
target aperture flux for each star. The final flux in instrumental units was then corrected to 
an infinite aperture (MIPS Data Handbook, V3.1, hereafter MDH3.1). 

Random uncertainties in the background-subtracted estimates of on-source flux were 
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estimated from the standard deviation of the mean of the multiple measurements for the 
24/zm observations. These "internal" uncertainties represent the precision of our measure- 
ments independent of the absolute flux calibration of Spitzer. The "internal" uncertainties 
for the 70/im and 160/xm observations of HD 12039 were estimated by the rms pixel-to-pixel 
dispersion inside the background annulus measured on the mosaic image and scaled to the 
area of the target aperture. 

Conversion to physical flux density units was performed with a simple multiplication 
of the instrumental total flux by the MIPS calibration factor as published currently in the 
MDH3.1. The MIPS absolute calibration uncertainties for the 24, 70, and 160 /xm wavelength 
bands are currently listed as 10%, 20%, and 20%, respectively in the Spitzer Observer's 
Manual Version 4.6 (hereafter SOM V4.6), and in the MDH3.1 3 . However, the median of the 
ratio of the measured 24/zm flux densities to Kurucz-model predicted 24/zm flux densities for 
the ensemble of 33 FEPS stars is 0.98 with a standard deviation 0.05 (§4.2). Assuming that 
the stars do not have infrared excesses, this small dispersion indicates that the MIPS 24/im 
absolute calibration is actually accurate to ~ 5%. This uncertainty estimate is consistent 
with 24/zm measurements of 69 F5-K5 (IV & V) stars observed by the MIPS GTO science 
team (Beichman et al. 2005a; Bryden et al. 2005). Excluding the sole warm debris system 
star HD 69830, Bryden et al. find that the mean MIPS 24yum to predicted 24/zm flux density 
ratio for the remaining 68 stars is 0.99 with a dispersion of 0.06, in excellent agreement with 
our results. We conclude that the photometric calibration of our MIPS 24/xm photometry 
relative to the best fit Kurucz models is 5-6%. Since we do not detect the photospheres 
at 70 and 160yum, we cannot independently derive uncertainties relative to Kurucz models. 
Therefore we adopt the SSC-published 20% absolute uncertainty for both the MIPS 70 and 
160/im photometry. Finally, we calculate a total uncertainty in the physical flux density 
for each photometric measurement by adding the "internal" and calibration uncertainties in 
quadrature. 

Low-resolution (R = 70-120) spectra were obtained on UT 2004 July 13 with the Infrared 
Spectrograph (IRS: Houck et al. 2004). We used an IRS high-accuracy, blue peak-up (la 
uncertainty radius = 0.4") to acquire the star in the spectrograph slit, thus minimizing slit 
losses and assuring high photometric accuracy. Two nod positions per cycle were obtained 
in standard staring mode with one cycle for the Short-Low wavelength range (7.4-14.5 /im), 
and 3 cycles for each of the two Long-Low wavelength ranges (14.0 - 21.3yum and 19.5 - 
38.0/im). The integration times were 6s per exposure. The spectrum beyond ~ 35/im suffers 
from high noise (Houck et al. 2004) and has been omitted. 



3 The MDH3.1 states that the absolute calibration at 24^m is "better than 10%." 
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The intermediate droopres products of the SSC pipeline Sll.0.2 were processed within 
the SMART software package (Higdon et al. 2004). The background was subtracted us- 
ing associated pairs of imaged spectra from the two nodded positions along the slit. This 
also subtracts stray light contamination from the peak-up apertures, and adjusts pixels with 
anomalous dark current relative to the reference dark frames. Pixels flagged by the SSC 
pipeline as "bad" were replaced with a value interpolated from an 8 pixel perimeter sur- 
rounding the bad pixel. 

The spatially unresolved spectra were extracted using a 6.2 pixel fixed-width aperture in 
the spatial dimension for the first order of the short wavelength low-resolution module, and 
5.1 and 3.1 pixels for the first and second order of the long wavelength low-resolution module, 
respectively. The spectra were calibrated using a spectral response function derived from IRS 
spectra and Kurucz stellar models for a set of 16 stars observed within the FEPS program 
that exhibit: 1) high signal-to-noise observations, 2) no residual instrumental artifacts, and 
3) no signs of infrared excess. The absolute flux density scale was tied to calibrator stars 
observed by the IRS instrument team, but reduced as for the rest of the FEPS sample and 
referenced to calibrated stellar models provided by the SSC (see also Hines et al. 2004a, 2005; 
Bouwman et al. 2005). Using this internally consistent approach enables us to improve the 
absolute flux calibration and reduces the noise on the adopted relative spectral response 
functions. This ensures that the uncertainties in the final calibration are dominated by 
photon noise and not by the uncertainties in the calibration. We estimate that the relative 
flux calibration across the spectrum is ~ 1 — 2%. As for MIPS, a direct comparison of the 
IRS spectra for the 33 stars to the predicted flux densities from our best fit Kurucz models 
indicates absolute uncertainties of « 6% at 24/im (§3). 

Infrared Array Camera (IRAC: Fazio et al. 2004) observations in channels 1, 2, and 4 
(3.6, 4.5, and 8.0/xm) were obtained on UT 2004 July 19 using the 32x32 pixel sub-array 
mode and a 4-point random dither pattern, with an effective integration time of 0.01s per 
image (frame time = 0.02s). The 64 images obtained at each of the four dither positions 
provided a total of 256 images of each star for a total integration time of 2.56 seconds in 
each channel. 

All photometry was performed on the Basic Calibrated Data products from the SSC 
S10.5.0 data pipeline as described in the SOM4.6 and the Pipeline Description Document 
available through the SSC. Aperture photometry was performed using IDP3 with a 3'/6 radius 
aperture centered on the target, and the background was estimated by the median of the 
pixels in a 12'.'2 - 39" radius annulus centered on the source in each sub-array image. This 
annulus circumscribes the 32x32 array and thus uses all pixels outside of a 10-pixel radius 
from the target. The background flux was normalized to the area of the target aperture and 
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subtracted from the summed target aperture flux. The reported source flux is the mean of the 
256 measures, corrected from a 3'/6 radius to the calibration aperture using values published 
in Table 5.7 of the Infrared Array Camera Data Handbook VI. 0. The "internal" uncertainties 
are estimated from the standard deviation of the 256 measurements, and are typically 
The ratios of measured to predicted flux densities for the 33 stars in the three IRAC bands 
are 1.08±0.04, 1.04±0.05, and 0.94±0.04 for IRAC 3.6, 4.5 and 8.0/im, respectively. This 
suggests that there are residual systematic offsets in the absolute calibration of the IRAC 
sub-array mode. These are probably caused by gradients in the filter band-passes across the 
detector (Quijada et al. 2004; Reach et al. 2005), which have yet been fully characterized 
for the sub-array mode. Given these systematic offsets, we conservatively adopt the 10% 
absolute uncertainties as listed for the three bands in SOM 4.6. 

The IRAC and MIPS photometry for HD 12039 is presented in Table 1. We also 
present photometry centered at 13, 24 and 33/im constructed from the IRS spectrum using 
rectangular band-passes of 1.6, 4.7 and 5/im FWHM respectively. These IRS flux densities 
are the error-weighted means with associated errors in the mean as appropriate. 

4. The Infrared Excess Associated with HD 12039 

Figure 1 shows a plot of the Fy(24/im)/F^(Ks) vs. F I/ (4.5/xm)/F ;/ (Ks) flux density ratios 
for the 33 stars in the sample. Data are shown from both the MIPS 24/im band and the 
synthetic IRS 24/im band. To construct the diagram, we used our most accurately calibrated 
IRAC band and the 2MASS i^s-band, which is least affected by extinction. In addition, since 
only three of the stars lie outside of the dust-free Local Bubble (d > 80 pc), and all are within 
150 pc, we do not expect extinction in the Ks-b&nd to cause us to identify spurious 24/im 
excesses. These ratios are independent of model atmosphere fits to the individual stars, 
and will identify excess candidates regardless of any relative calibration offsets between the 
instruments. 

Figure 2 shows the MIPS F„(24/im) /F 1/ (Ks) ratio plotted versus the IRS F J/ (24/im) /F„(Ks) 
ratio illustrating the independent instrumental calibration for IRS and MIPS as compared 
with 2MASS. The 30 Myr old, G3/5V star HD 12039 4 is clearly distinguished in Figures 1 
& 2, and apparently has significant excess emission at 24/im detected in both IRS and MIPS 
independently relative to the other 32 stars. 

Since the MIPS and IRS data for each object were obtained ~ 2 days apart, the agree- 



Dctails concerning our age determination for this star are given in Appendix A. 



- 8- 



ment between the flux density measurements for HD 12039 rules out a chance superposition 
of an asteroid in one of the observations as the explanation for the excess. Furthermore, 
inspection of 2MASS images shows that there are no other objects having Ks < 15 mag 
within 5" of the position of this star, and the 24/xm image is indistinguishable from a single, 
isolated point source. Therefore the measured 24/im emission is associated directly with 
HD 12039. The excess in the F J/ (24/im)/F 1/ (Ks) ratios are observed at the 4.2 and 4.3a level 
relative to the median ratios for the other 32 stars from the MIPS 24/im and IRS 24/im 
photometry (Figs. 1 & 2). The probability of drawing this > 4<r result from a Gaussian 
distribution of data for these 33 targets is < 33 x 1.3 x 10 -4 = 4 x 10~ 3 . 

As a further check on the reality of this excess at 24/im, we compared the distributions of 
of the ratios our measured MIPS and IRS 24/im flux densities to the flux densities predicted 
by Kurucz model atmospheres fit to optical and near-IR photometry for all 33 stars (Figure 
3). Excluding HD 12039, the median and standard deviations (la) of these distributions are 
0.98 ± 0.05 and 0.99 ± 0.06 for MIPS 24/im and IRS 24/im respectively, suggesting that 
the MIPS and IRS absolute flux calibrations are accurate to 5-6% relative to the best fit 
Kurucz models. The ratio of measured-to-Kurucz model 24/im flux densities obtained with 
both instruments for HD 12309 is 1.33±0.04, and [F„ j0bs (24/im) - F !/imodel (24/im)]/a is 5.4 
and 5.3 for MIPS and IRS, respectively, where a is the rss of the standard deviation of 
the distributions in Figure 3 and the total measurement uncertainty for HD 12039. Based 
upon the above analysis, we conclude that HD 12039 has excess emission above its stellar 
photosphere at 24/im. 

Figure 4 shows the IRAC and MIPS photometry plus the IRS low resolution spectrum 
of HD 12039 compared to the model photosphere. The object was not detected at 70/iin 
with a 3o" upper limit that is 6 times the expected photospheric flux. Upper limits are also 
shown for 160/im, 1.2mm, and 3.1mm. 

The photospheric emission component of the spectrum in Figure 4 was fit with Ku- 
rucz model atmospheres (including convective overshoot) to the optical and near-infrared 
photometry. Predicted magnitudes were computed as in Cohen et al. (2003 and references 
therein) using the combined system response of detector, filter and atmosphere (for ground- 
based observations). The best-fit model was computed in a least squares sense with the 
effective temperature and angular diameter of the star as free parameters, [Fe/H] fixed to 
solar metallicity, surface gravity fixed to the value appropriate for the adopted stellar age 
and mass (log g = 4 in this case), and visual extinction fixed to mag within the dust-free 
Local Bubble (d* ~ 42 pc). The derived best-fit stellar model has T e ff = 5688K, consistent 
with HD 12039's catalog G3/5V spectral type. 

In comparing Spitzer data to the model photosphere we have adopted the weighted 
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average wavelengths that correspond to the monochromatic flux density within the filter 
band-passes according to the SOM4.6. Also, for the purposes of plotting the data, a color 
correction (F c = 1.056 x F ^, s ) has been applied to the MIPS 24/zm point assuming a color 
temperature T = 110K estimated from the slope of the IRS spectrum from 24/xm to 33/im 
(see below) 5 . 

The IRAC and short wavelength IRS data for HD 12039 are consistent with the best-fit 
Kurucz model. However, the IRS spectrum begins to depart from the photosphere between 
12-14/im, passes through the MIPS 24/xm datum, and continues to increase until at least 
33/im. This is shown more clearly in the bottom panel of Figure 4, where we have plotted 
the ratio of the measured flux densities to the model-predicted flux densities. We also show 
the synthetic 33/tm photometric point from the IRS spectrum, which lies > 3.5a above the 
Kurucz model. 

The lack of significant excess continuum emission for wavelengths < 14/im and the 
measured 24/im excess together imply a maximum color temperature of T c < 160K. The 
slope and uncertainty between the MIPS 24/im and IRS 33/im broadband photometric points 
correspond to color temperatures T c = IO9+37K, ~ 110K. The 70/im upper limit constrains 
the 33 — 70/im color temperature to > 80K. 

In the next section we further characterize the infrared excess, and explore circumstellar 
disk models. We also consider constraints on the nature of the thermally emitting material 
imposed by timescales for evolution of dust in the system. 

5. Analysis 

5.1. General Characteristics of the Dusty Debris System 

To further characterize the infrared excess, we assume that the emission originates from 
large grains that absorb and emit radiation efficiently at all relevant wavelengths; that is, 
we treat the grains as blackbodies (e.g., Aumann et al. 1984; Walker & Wolstencroft 1988; 
Backman & Paresce 1993). Blackbody emission at wavelengths as long as 40/tm requires 
grain radii larger than a ~ 7/tm, and for such grains the dust physical temperature T dust is 
the same as the color temperature. 

Given the fitted color temperature of T ~ 110K and assuming blackbody grains in 
thermal equilibrium with the stellar radiation field, we calculate a typical grain location of 6 



5 MIPS color-corrections are tabulated in the MIPS Data Handbook V3.0 available from the SSC. 
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AU from the star. We also derive an upper limit to the grain radial distance by considering 
the 3cr upper limit to the 70/im flux density (Figs. 3 & 4). Only blackbody grain models with 
R out < 11AU do not violate this constraint. Thus, the observed infrared excess spectrum for 
HD 12039 is consistent with a fairly narrow range of dust temperatures and radial annuli. 
We note that although the IR emission from HD 12039 may be dominated by such relatively 
large grains, the characteristic emission temperature of 110K is low enough that the lack of 
solid state emission features at e.g. 10 and 18/xm from the bending and stretching modes of 
standard astronomical silicates does not exclude the presence of smaller grains. 

The observed excess at 24/xm and limit at ~ 14/xm excludes the presence of material at 
T > 200K with more than ~ O.lx the surface density inferred for the material near 110K, 
assuming that any material at higher temperatures and smaller radii from the star would 
have a surface density oc r° as for Poynting-Robertson (P-R) drag. Similarly, the 70/zm limit 
excludes material at T~ 50K with more than ~ 0.5 x the total grain cross-sectional area of 
material at 110K. 6 The fractional infrared luminosity of the excess across the wavelength 
range 24-70/im is f = Ltr/L* = 1 x 10~ 4 . 

5.2. Disk Models 

We further investigate the HD 12039 debris system using the dust disk models of Wolf & 
Hillenbrand (2003), which take into account absorption and emission from dust grains based 
on their optical properties as opposed to the simple blackbody assumption. In particular, 
we use "astronomical silicates" (Draine & Lee 1984; Laor & Draine 1993; Weingartner & 
Draine 2001). We assume a volume density profile n(r) oc r _1 that corresponds to a disk 
with constant mass surface density S(r) oc r°, consistent with but not requiring dynamic 
control by P-R drag (§5.3). 

As discussed above, the wavelength at which the dust re-emission departs significantly 
from the stellar photosphere (A ~ 14 /jm) yields an initial estimate of the inner radius of 
the disk R in . This is also the location of the warmest material, but the dust temperature 
of T~ 110K is too low to produce detectable silicate emission features (see, e.g., Fig. 10 in 
Wolf & Hillenbrand (2003) and their discussion of the effects of inner holes, thus maximum 
temperatures, on the predicted spectra of debris disks). 

To explain the spectrum of the infrared excess in HD 12039, we calculated the inner and 



6 In this case we can only place limits on dust emitting at a single temperature; there are no constraints 
on the surface density. 
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outer disk radii and corresponding total dust mass for a suite of single-size-grain models with 
particle radii a = 0.1,0.5,1,3,7,25,50 and 100/rni. We applied the Levenberg-Marquardt 
algorithm (Marquardt 1963; Markwardt 2003) to find the models that best fit the observa- 
tions of HD 12039 in a chi-squared sense (Rodmann et al. in prep). The parameters of these 
models are presented in Table 2 and illustrated graphically in Figure 5. Models assuming 
small grains require large orbital radii from the star (i? in R3 10-40 AU, R out ~ 20-50 AU), 
because these small grains have low radiative efficiencies compared with blackbody grains. 
For grain sizes a > 7//m, we find values of R in = 4 AU and R out = 6 AU, approximately in- 
dependent of particle size, consistent with the simple blackbody dust model. Table 2 and 
Figure 5 also show that a model assuming a = 7/zm grains (located between 4-6AU) yields 
the minimum dust mass. These single-size-grain models illustrate the inherent degeneracies 
in fitting spectral data alone without independent morphological information such as direct 
imaging of the systems in thermal emission or scattered light. 

Real systems in which the dust is produced by collisions of parent bodies will likely 
be composed of a distribution of particle sizes (Dohnanyi 1969; Williams & Wetherill 1994; 
Tanaka et al. 1996). Therefore, we also modeled the infrared excess for HD 12039 us- 
ing a power-law grain size distribution n(a) oc a~ p that included small grains of radius 
a min « 0.4 /im (just smaller than the blowout size, see below) and a maximum grain size 
of a max = 1000 /zm (just larger than can efficiently radiate at Spitzer wavelengths). The 
power-law exponent was set to p = 3.5, which is produced by a collisional cascade. For such 
a distribution, most of the opacity is supplied by the smallest grains and most of the mass 
in is contained in the largest grains. For these assumptions the best-fit inner and outer disk 
radii are R m = 28 AU and i? ou t = 40 AU, respectively (Table 2). In general, more mass is 
required by this model than by the single-size-grain models considered above because of the 
inclusion of very large grains with small ratios of emitting area to mass. 

Disk models and the photosphere-subtracted Spitzer data for HD 12039 are shown in 
Figure 6. All of the models require an inner hole radius R m > 4AU, but we are unable to 
choose between the various models from the spectral fitting alone. We therefore consider, 
in the next section, physical scenarios for dust production and migration that lead us to a 
preferred model. 

5.3. Estimating Time Scales for Grain Evolution 

We further constrain properties of the dust debris system around HD 12039 by con- 
sidering the interaction between dust grains and radiation from the star (e.g., Burns et al. 
1979; Backman & Paresce 1993; Wyatt et al. 1999, and references therein; Thebault et al. 
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2003). Assuming there is no remnant gas in the system, typical astrophysical silicates with 
radius a < 0.5pm and density p = 2.5 g cm -3 produced by parent bodies in circular orbits at 
any distance surrounding a star with the luminosity and mass of HD 12039 will be expelled 
by radiation pressure. When grains travel a factor of four farther from the star than from 
where they originated, their equilibrium temperature falls by a factor of two and they are 
no longer important contributors to the observed infrared emission. The relevant timescales 
for this movement are ~ 8yrs and ~ 80yrs for grains starting at 6AU and 30AU respectively. 
Ice grains (p ~ 1.0 g cm" 3 ), assumed to be "dirty" with albedos similar to silicates, will be 
similarly removed if a < 1.5pm. Therefore, the presence of grains smaller than the blowout 
size would signal a very recent collisional event. 

In addition to radiation pressure expelling small grains from the system, P-R drag 
will cause inward migration of particles in the absence of remnant circumstellar gas. The 
time scales for silicate grains (p = 2.5 g cm~ 3 ) to spiral into the star from a release radius 
^dust ~ 6AU are tp-p ~ 7 x 10 4 or ~ 5 x 10 5 yr for grains with a ~ 1pm (a bit above the 
blow-out size) or a ~ 7pm (the minimum blackbody emitter size), respectively. This is much 
longer than the blow-out time (< 10 2 years), but still short compared to the t ~ 30 Myr age 
of the system. Similar sized grains at 30-40AU would have P-R timescales ~ 1 — 10 Myr, 
again shorter than the system age. In either distance regime, small particles are removed 
on sufficiently short timescales that the material must be replenished if we are witnessing 
a steady-state. Alternatively, as mentioned above, the dust could have been produced by a 
relatively recent event. 

The dust fractional surface density (m 2 cross-section per m 2 of disk) near the minimum 
orbital radius R m i n is numerically roughly equal to the fractional luminosity in the excess f = 
LipjL^ (Backman 2004), i.e. f~ 1 x 10~ 4 . From this we can explore the relative importance 
of P-R drag compared with grain evolution caused by mutual collisions (e.g., Wyatt et al. 
1999; Wyatt & Dent 2002; Wyatt 2005). Considering our two primary disk models discussed 
above, the dust mutual collision timescales in this system are approximately 2 x 10 4 yr for 
7pm grains at r^st = 6AU, and ~ 2 x 10 5 yr for 0.5pm grains at rd us t = 30AU with scaling as 
r 1 ' 5 (Backman & Paresce 1993). At 30AU, the collision and P-R timescales are comparable, 
so material moves towards the star from its region of origin before collisions can modify it 
significantly. At 6AU the collision timescale is 25 times smaller than the P-R timescale, 
so collisions can grind material down to the blowout size before the larger grains migrate 
towards the star via P-R drag. This creates a natural paucity of material inside the inner 
radius of the parent body orbits, and could explain the observed lack of significant material 
at hotter temperatures (i.e., smaller orbital radii). This mechanism for keeping inner holes 
clear may operate in other well studied disks (e.g., Wilner et al. 2002; Dominik & Decin 
2003; Wyatt 2005). 
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In this analysis we have not considered corpuscular drag from collisions with the stellar 
wind from HD 12039, which would increase the minimum grain blowout size. Even for young 
Sun-like stars, corpuscular drag forces are thought to be small compared with radiation forces 
(Fwind/Fp_R ~ 0.35: Gustafson 1994). Very tiny grains (e.g., a < 0.05/zm for silicates) would 
be inefficient absorbers at the peak wavelengths of the stellar spectrum and perhaps stable 
against radiation blowout, but possibly vulnerable to ionization and resulting magnetic forces 
in ways that are poorly understood. We can assume that silicate grains with radii at least 
over the range 0.05 < a < 0.5/iin (or dirty ice grains with 0.1 < a < 1.5/im) are unstable 
against ejection from this system. 

6. Discussion 

6.1. The Debris System 

Our inferred properties of the debris associated with HD 12039 assume that dust dy- 
namics are dominated by interactions with the stellar radiation field and are not determined 
by dust-gas dynamics. This requires the gas-to-dust ratio to be < 0.1 (e.g., Takeuchi & 
Artymowicz 2001; Klahr & Lin 2001). While this is consistent with the known properties of 
HD 12039 (no signatures of active accretion in the stellar spectrum, no strong mid-IR atomic 
or molecular emission-line features), detailed searches for remnant gas have not been made. 
With this caveat, the short grain evolution timescales relative to the system age imply that 
the dusty debris consists primarily of second-generation grains released by collisions and/or 
sublimation of larger parent bodies. 

We have presented a range of possible models that successfully fit the thermal excess, 
but two are of particular interest because they represent end points of a plausible range of 
grain physical characteristics: (1) large "blackbody" grains (a > 7/im) located between 4 
and 6 AU from the star, and (2) grains with sizes 0.4 < a < lOOOyum characterized by a 
powerlaw distribution, and located between 28 and 40 AU. Without resolved images of the 
debris system or more sensitive searches for mineralogical features, we cannot distinguish 
between these models with confidence. However, the surface area in model (2) with the wide 
belt and a grain size distribution that is dominated by grains near the blowout size would 
require either: a) a large mass (> 100 Msarth) of parent bodies in collisional equilibrium to 
maintain the small-grain population against rapid blowout and P-R drift; or b) the discovery 
of a large transient dust signature observable only for a few hundred years and due to recent 
collisional event. 

Based on planetesimal accretion code models of a minimum-mass solar nebula (MMSN) 
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surrounding a Sun-like star, Kenyon & Bromley (2004, 2005) have shown that the growth 
of planetesimals in the terrestrial zone can produce ~ 1 and ~ 3 mag infrared excesses at 
10/im and 24/im respectively. The peak emission occurs from ~ 10 4 to ~ 10 6 yrs when the 
largest planetesimals have reached ~ 2000 km in size. By 10 7 yrs, the excess is nearly gone 
at 10/im, but persists at 24/im until ~ 10 8 yrs. Importantly, even though the overall dust 
production (thus IR emission) declines, individual collisions between 10-100 km bodies will 
produce short-lived ~ 0.5 — 1.5 mag signatures (see Fig. 3 in Kenyon & Bromley 2004; 
Kenyon & Bromley 2005). The observed ~ 0.3 mag excess at 24/im in HD 12039 is roughly 
consistent with this scenario. In contrast, at it is harder to grow large bodies at farther 
distances from the star. Collisions at 30 AU are predicted not to be massive enough in 
general to produce an observable excess above the photosphere at 24/im (Kenyon & Bromley 
2005). These considerations suggest that our model (1) is the most plausible explanation for 
the infrared excess in HD 12039. 

The lack of material much warmer than ~ 110K in the HD 12039 system could be 
explained by several mechanisms. First, this temperature is approximately equal to the 
threshold for the onset of rapid sublimation of micron-sized water ice grains. Thus, if the 
dust parent bodies are primarily icy, the inner edge of the observed debris distribution could 
be determined by grain evaporation. Second, as discussed in §5.3, the timescale on which 
grains collide and produce fragments small enough for rapid blowout is an order of magnitude 
shorter than the P-R timescale at 6 AU, but the timescales are comparable at 30 AU. This 
suggests that the inner belt model could be self-limited. Dust grains would collide 'in situ' 
before drifting inward a significant distance from their production location (Wyatt 2005), and 
the small grains would be blown out by radiation pressure. The existence of a planetesimal 
belt surrounding HD 12039 might require the presence of a nearby giant planet to stop the 
planet formation process, as is understood to be the case with Jupiter and the asteroid belt 
in our solar system. Third, a planet could exist at the inner edge of the debris annulus and 
consume/scatter inbound grains (e.g., Moro- Martin, Wolf & Malhotra 2005). 

At the distance of HD 12039, and with Lpr/L* « 1 x 10~ 4 in the range of objects 
detected previously in scattered light with HST (e.g., Schneider et al. 1999; Ardila et al. 
2004; Fraquelli et al. 2004), the debris disk should be easily resolved in scattered light by 
coronagraphic imaging with NICMOS aboard HST, if it is composed of small grains at 30AU. 
A null detection would strengthen the case that the material is concentrated in a belt at 
4-6AU. NICMOS (and ground-based AO) imaging should also place limits on the presence 
of a planet perturber that may be responsible for stirring up the debris. 
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6.2. Terrestrial Planetary Debris Systems 

HD 12039 is the only object among the 33 Sun-like stars in this FEPS sub-sample that 
exhibits a prominent debris system at relatively warm temperatures approaching those of 
asteroidal material in our own solar system. It also apparently lacks significant material at 
cooler, Kuiper Belt-like temperatures. We know of only a handful of other such "Terrestrial 
Planetary Debris Systems" (Table 3). 

Three systems are associated with stars more massive, and hence more luminous, than 
the Sun. The 8 Myr old AOV star HR 4796A shows a strong infrared excess that is well 
modeled by a single temperature blackbody (T = 110K) with grains orbiting at ~ 35AU 
(Jura et al 1998). Further high resolution imaging of HR4796A in the infrared indicates the 
presence of a warmer T = 260 ± 40 K, compact dust component located at 4 AU (Wahhaj et 
al. 2005). The 10-20 Myr old F3V star HD 113766A is part of a wide (170 AU) binary that 
is a kinematic member of Lower Centaurus Crux association. The T~ 350K material has 
been modeled by a debris belt 0.4-6 AU from the star (Meyer et al. 2001; Chen et al. 2005). 
The 300 Myr old A2V star ( Lep was first identified by Aumann & Probst (1991) to possess 
an infrared excess, and Chen & Jura (2001) find dust of several hundred K, corresponding 
to a spatial distribution <6-9 AU. 

Two additional warm debris systems are associated with stars less massive than the Sun: 
HD 98800 and HD 69830. The 8-10 Myr old quadruple system HD 98800 in the TW Hya 
association is comprised of two spectroscopic binaries with K and M primaries (Torres et al. 
1995; Soderblom et al. 1996, 1998). A circumbinary debris system first detected by IRAS 
(Walker & Wolstencroft 1988; Zuckerman & Becklin 1993) orbits HD 98800B (Koener et al. 
2000) is well characterized by a single T^ 160K blackbody spectrum from 8/im to 7mm, 
implying debris orbiting at 2-4AU (Low et al. 1999; Hines et al. 2004b; Low et al. 2005). 
The presence of silicate emission features suggest that there is an additional population of 
small grains in this system (Skinner et al. 1992; Sitko et al. 2000; Schiitz et al. 2005). 

The 2 Gyr old K0V star HD 69830 was also identified by IRAS as harboring a debris disk 
system (Mannings & Barlow 1998). Similar to HD 12039, Spitzer observations of HD 69830 
reveal a strong 24/xm excess, and the 70/im flux density is consistent with the photosphere 
(Beichman et al. 2005b). In contrast to HD 12039, however, the IRS spectrum of HD 69830 
has abundant mineralogical features similar to those found in ISO and ground-based mid-IR 
spectra of comet Hale-Bopp (Crovisier et al. 1996, 1997; Wooden et al. 1999a,b). Beichman 
et al. (2005b) model the system with a = 0.25/iin grains at about 3AU, and suggest a recent 
dust-producing event. 

Finally, the 30Myr old G3/5 star HD 12039 and 300Myr old GOV star HIP 8920 (Song 
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et a. 2005) each have masses and luminosities very similar to the Sun. HIP 8920 is extremely 
unusual, exhibiting a very hot (T^ ust = 650K) bright (Ltr/L* ~ 0.04) infrared excess. As for 
HD 69830, but unlike HD 12039, HIP 8920 shows a strong lOyum silicate emission feature 
signaling the presence of grains smaller than 3/im. Song et al. (2005) note these small, hot 
grains would reside at ~ 0.4AU. 

Twenty isolated Sun-like stars in the 100 Myr old Pleiades have been observed by the 
FEPS program to look for excess infrared emission (Stauffer et al. 2005). One "probable" 
and two "possible" detections at 24/im were found. These qualifiers result from the fact the 
Pleiades members are passing through a molecular cloud, and are embedded in ISM dust 
that can be illuminated by the stars. Thus they may have an associated infrared excess that 
is unrelated to the debris disk phenomenon. Also, the IRS 33/zm and 70/zm upper limits 
suggest that these disks, if real, are fairly cool T< 85K (Stauffer et al. 2005). Therefore, we 
have not explicitly included these stars in our list of warm debris disks. 

The frequency of warm debris systems associated with Sun-like stars may be very 
coarsely estimated based on Spitzer results presented to date. So far, FEPS has found 
1/33, Chen et al. (2005) find 1/40 (but possibly 3 total), Beichman et al. (2005a, 2005b) 
and Bryden et al. (2005) have found 1/84. A MIPS GTO investigation of the lOOMyr old 
open cluster M47 found 2/17 Sun-like stars with a 24/zm excess, one of which is sufficiently 
strong that it might be from terrestrial temperature dust (Gorlova et al. 2004). Another 
Spitzer investigation of 24 of the known stars within the TW Hydra Association (TWA: 
Low et al. 2005) confirmed that HR 4796A and HD 98800B are the only members of the 
association to harbor warm debris disks 7 . As discussed above, FEPS observations of the 
Pleiades may have also uncovered 3/20 objects with 24/zm excess. 

Compiling these results suggests that approximately ~ 1 — 5% of Sun-like stars in the 
age range ~ 8 Myr - 300 Myr (plus HD 69830 at 2 Gyr) exhibit infrared excesses peaking 
at terrestrial temperatures. This rarity may reflect the short duty cycle of such events or, in 
steady state, an intrinsically low frequency of stars with the observed level of dust production 
in terrestrial planet-building zones. Kenyon & Bromley (2005) find that the rate of collisions 
sufficiently large to produce a substantial 24/xm signal is < 10~ 5 yr _1 for the period of 
terrestrial planet formation t ~ 1 — lOOMyrs 8 . Therefore the apparent paucity of Sun-like 



7 A tentative (2a) detection at 12 and 18^m of a T=170K disk in TWA 17 by Weinberger et al. (2004) 
was not confirmed by MIPS 24^m observations that are consistent with the stellar photosphere (Low et al. 
2005). 

8 In our own solar system, the terrestrial planets are thought to have been 80% complete by an age of 30 
Myr (Kleine et al 2002, 2003; Yin et al. 2002). 
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stars with bright thermal excesses and terrestrial temperatures may not be surprising. 

7. Conclusion 

We have discovered a debris system around the Sun-like star HD 12039 that exhibits 
some properties analogous to the dust associated with the asteroid belt in our solar system. 
While very small particles (a < 1/im) at Kuiper Belt distances (r = 30 AU) cannot be ruled 
out completely, we find that a debris belt consisting of a > 7fxm grains orbiting HD 12039 
at rdust = 4 — 6AU provides the most plausible geometry. As for our solar system, such a 
narrowly confined zone of material may signal the presence of (as yet unseen) planets in the 
system. 

HD 12039 joins only a handful of other debris disk systems with similar properties 
out of hundreds observed, suggesting that such systems are rare. If this rarity reflects the 
true incidence of terrestrial-zone debris, then perhaps asteroid belt analogs are uncommon. 
Alternatively, the collisional events that produce such strong infrared signals may be rare. 
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Center/ California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. FEPS is pleased to acknowledge sup- 
port from NASA contracts 1224768, 1224634, and 1224566 administered through JPL. The 
MPIA team is supported by the EU-RTN-HPRN-CT-2002-00308 'PLANETS' network and 
the German Research Foundation (DFG) through the Emmy Noether grant WO 857/2-1. 

A. The Age of HD 12039 

Song, Zuckermann & Bessell (2003) claim that HD 12039 (DK Cet, HIP 9141) is a mem- 
ber of the ~30-Myr-old Tucana-Horologium association (d* = 42±2pc, Perryman 1997). 
We independently confirmed the kinematic association by comparing the proper motion 
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(Zacharias et al. 2004) and mean radial velocity (5.0 ±0.2 km s" 1 ; Song et al. 2003; Nord- 
strom et al. 2004; White et al. 2006) to the values expected for a star at d = 42 pc with the 
Tuc-Hor group space motion (Zuckerman & Webb 2000). 

An age of ~30 Myr is also consistent with other estimates. Position in the HR diagram 
(T e// =5688±60K, log(L/L ) = -0.05±0.04 dex: Carpenter et al. 2005) compared to evo- 
lutionary tracks suggests an age of 28 Myr (D'Antona et al. 1997) to 33 Myr (Baraffe et al. 
1998); the mass predicted from either set of tracks is l.O2M . 

Chromospheric activity and lithium depletion data corroborate the isochronal age es- 
timate. The Call emission index has been measured to be log R HK = -4.21 (Soderblom, 
private communication) and log R#x = -4.14 (White et al. 2006). This is broadly consistent 
with an age of 30 Myr although the uncertainties in estimating ages from Call indices for 
ages < 300 Myr are large. The Li I A6707 equivalent width (Song et al. 2003; Wichmann 
et al. 2003; White et al. 2006) when interpolated among the distribution of Li I values for 
Sun-like cluster members (Mamajek, in preparation) suggests an age of 15^^ Myr. Consid- 
ering all of these age indicators, we adopt an age of ~30 Myr for HD 12039 based largely on 
its kinematic membership in the Tuc-Hor association. 
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Fig. 1. — Flux ratio diagram [F„ (24/mi) /F„ (Ks) vs. F I ,(4.5^m)/F„(Ks)] for the 33 stars in the FEPS 
VI. 1 data release. Data are shown from both the MIPS 24/im broad band photometry and from a syn- 
thetic broad band point centered at 24/im formed from the IRS spectra. Note that the F l/ (4.5/^m)/F, y (Ks) 
ratio for HD 12039 is near the median for the other 32 stars consistent with photosphcric colors, but the 
F„(24/im)/F„(Ks) ratio is anomalously high indicating a significant 24/j.m excess above the photosphere. 
Error bars represent la in the ratios calculated from the internal IRS and MIPS measurement uncertainties, 
and the total uncertainties of the 2MASS K s flux densities. The 24/im excess from HD 12039 is detected in 
both MIPS and IRS independently and at different epochs. The dotted lines indicate the ratios predicted by 
the best- fit Kurucz model atmosphere for HD 12039. Stars with photosphcric temperatures different than 
HD 12039 will have slightly different F„(4.5/im)/Fi,(Ks) ratios, which explains the dispersion in this quantity 
for the 33 stars in the sample. 
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Fig. 2.— The F^(24/im)/F^(Ks) ratio for measurements from IRS & MIPS for the 33 stars 
in the FEPS VI. 1 data release. HD 12039 stands out clearly as having an anomalously high 
F l/ (24/xm)/F I/ (Ks) ratio in both instruments indicating a significant 24/zm excess above the 
photosphere. The solid line is of slope = 1. 
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Fig. 3. — Histograms of the ratios of the measured IRS (top) and MIPS (bottom) 24/im 
flux densities to the predicted flux densities from our best-fit Kurucz models for the 33 stars 
in the FEPS VI. 1 data release. The bin width represents the average uncertainty in the 
ratios, and the crosses above the histograms indicate the median and standard deviation of 
the measurements excluding HD 12039 (1.00 ± 0.06 and 0.98 ± 0.05 for IRS and MIPS, 
respectively). As in Figures 1 & 2, HD 12039 stands out clearly as having a significant 24/im 
excess above the photosphere, and is detected in both MIPS and IRS independently and at 
different epochs. 
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Fig. 4. — Spectrum of HD 12039. The upper panel shows the ground and space-based 
photometry (solid symbols), the IRS spectrum (green curve), and the Kurucz model that 
best fits the optical and near-IR photometry (solid black curve). The MIPS 24/im point has 
been color-corrected assuming a color temperature T = 110K (see text). The la error bars 
for each point are plotted, but they are often smaller than the dot size in the top panel. 
We also show upper limits from MIPS 70 and 160/im photometry, as well as upper limits 
at millimeter wavelengths (Carpenter et al. 2005). The upper limits represent the actual 
on-source measured flux densities minus the background (which can be negative) plus three 
times the total uncertainty including the absolute calibration uncertainty. Also shown is the 
best-fit emission model for simple blackbody grains (dot-dashed line), and a blackbody model 
with the outer radius extending to 11AU (dotted line). Note that this latter model violates 
the 3cr upper limit at 70/zm and thus is ruled out by the observations. The lower panel shows 
the spectrum of HD 12039 divided by the Kurucz model. The error bars for the IRS 33fim 
synthetic photometric point are duplicated to the right of the actual point for convenience. 
The IR-emission can be seen departing from the photosphere for A > 12 — 14/im. 
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Fig. 5. — The single- size- grain debris disk models for HD 12039 following the prescription 
outlined in Wolf & Hillenbrand (2003). The models are discussed in §5.2 and their parameters 
are listed in Table 2. The inner and outer radial distances of the disk from HD 12039 and 
the associated dust mass for each disk model are plotted as a function of grain size. The 
vertical dotted line labels the blowout size (a = 0.5/im) for silicate grains; smaller grains 
are expelled from the HD 12039 system by radiation pressure on timescales < few hundred 
years. Note that grains with a > 7/xm behave as blackbodies and orbit in a belt from 4-6AU. 
The model with a = 7/im grains represents the minimum dust mass. 
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Fig. 6. — The spectrum of HD 12039 compared with spectra produced by the disk models 
from Table 2, all divided by the best-fit Kurucz model atmosphere. Several single- size-grain 
models are depicted as well as a model composed of grains with a powerlaw distribution of 
sizes from 0.4-1000/im (Table 2; see text). The best fit blackbody model is also shown for 
comparison. The single-size-grain dust models with particle sizes a > 7/im from Table 2 are 
indistinguishable from the simple blackbody model, and are not plotted independently in 
this figure. 
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Table I. Spitzer Photometry of HD 12039 



Band 


Wavclcngth a 
H 


Flux Density 
(mJy) 


fint 

(mJy) 


(mJy) 


Flux Density Excess 
(mJy) 


Instrument 


3.6 


3.535 


777.50 


1.43 


31.10 


71.51 


IRAC 


4.5 


4.502 


456.17 


1.61 


22.81 


26.17 


IRAC 


8.0 


7.735 


158.72 


1.03 


6.35 


-7.28 


IRAC 


13 d 


13.17 


61.65 


0.16 


9.25 


4.3 


IRS 


24 


23.68 


24.65 


0.22 


1.50 


6.1 


MIPS 


24 d - e 


24.01 


24.63 


0.38 


1.53 


6.0 


IRS 


33 d 


32.36 


20.31 


0.23 


3.05 


10.5 


IRS 


70* 


71.42 


4.5 


2.5 


2.7 




MIPS 


160 f 


155.9 


-42 


23 


28 




MIPS 



a Weighted average wavelengths (SOM 4.6). 

b The total uncertainty (atot) is the rss of the "internal" measurement uncertainty (precision) and 
the "calibration" uncertainties (accuracy). 

c Calculated relative to the photospheric model described in the text. 

d Synthetic band-pass photometry with 1.6, 4.7 and 5^tm band-widths (FWHM) for 13, 24 and 33^tm 
respectively. 

°The MIPS 24fim photometry has not had a color-correction applied. The current value for the color- 
correction would yield F c = 1.056 X F a ^ s = 25.94 mJy, so the thermal excess above the photosphere 
would be 7.4 mJy). 

f Upper limits for the purposes of model fitting and plotting were calculated by F (upper) = 
F(observcd) + 3<7tot- 



Table 2: Model Parameters for the Debris System HD 12039 



Model 


Grain radius 
(H 


Inner disk 
radius (AU) 


Outer disk 
radius (AU) 


Total mass 
10- 4 M carth 


la 


0.1 


40 


50 


0.26 


lb 


0.5 


27 


40 


0.1 


lc 


1 


20 


31 


0.08 


Id 


3 


11 


19 


0.07 


le 


7 


4.0 


6.0 


0.02 


If 


25 


3.7 


6.0 


0.1 


lg 


50 


3.8 


6.1 


0.2 


lh 


100 


3.9 


6.0 


0.4 


2 


0.4-1000 


28 


40 


4.5 



Table 3. Terrestrial Planetary Debris Systems 



Star 


Sp. Type 


age 


Tdust 


r dust 


Lir/L* 


Reference 






(Myrs) 


(K) 


(AU) 


(xlO- 4 ) 




HR 4796A 


AOV 


8-10 


110 


35 


50 


1 


C Lcp 


A2V 


300 


370 


< 6 


1.7 


2 


HD 113766A 


F3V 


15 


350 


0.4-6 


21 


3 


HIP 8920 


GOV 


300 


650 


0.4-1 


400 


4 


HD 12039 


G3/5V 


30 


110 


4-6 


1 


5 


HD 98800B 


K7V 


8-10 


160 


2-4 


1900 


6,7 


HD 69830 


KOV 


2000 


207 


1 


3 


8 


Sun 


G2V 


4500 


150-170 


1-4 


0.001 


9,10,11 
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